Non-Markovian spin relaxation in two-dimensional electron gas. 



M.M. Glazov 

F. Ioffe Physico-Technical Institute, Russian Academy of Sciences, 194021 St. Petersburg, Russia 



o 
o 

(N 
Or 

o 



E.Ya. Sherman 

Department of Physics and Institute for Optical Sciences, University of Toronto, 
60 St. George Street, Toronto, Ontario, Canada M5S 1A7 

We analyze by Monte-Carlo simulations and analytically spin dynamics of two-dimensional elec- 
tron gas (2DEG) interacting with short-range scatterers in nonquantizing magnetic fields. It is 
shown that the spin dynamics is non-Markovian with the exponential spin relaxation followed by 
the oscillating tail due to the electrons residing on the closed trajectories. The tail relaxes on a long 
time scale due to an additional smooth random potential and inelastic processes. The developed 
analytical theory and Monte-Carlo simulations are in the quantitative agreement with each other. 
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Spin-dcpcndcnt phenomena in semiconductors and 
semiconductor nanostructurcs attract an increasing in- 
terest during last decade (for review see, e.g. Ref. QJ). 
The understanding of the mechanisms of spin decohcr- 
ence, usually considered as a Markovian process, as well 
as the possibilities to control it is an important prob- 
lem in the field. In the two-dimensional (2D) zinc-blende 
structures the electron spin dynamics is governed by the 
spin-orbit (SO) splitting of the conduction band. It orig- 
inates from the lack of the inversion center either in the 
bulk material as represented by the bulk inversion asym- 
metry (BIA or Dresselhaus term 0,0) or the structural 
asymmetry of the heteropotcntial revealed by the Rashba 
term . In all the cases the SO interaction gives rise to 
the effective Zeeman magnetic field which is characterized 
by the electron wavevector fc-dependent spin precession 
vector Clk whose direction determines the axis and the 
absolute value is the spin precession rate. In the colli- 
sion dominated regime of electron motion with fir <C 1 
(where is the typical value of \ flk\ and r is the scat- 
tering time) the small spin rotation angles f2r between 
successive collisions arc not correlated, thus the spin re- 
laxation rate t^ 1 is proportional to the f2 2 r. 0,01 This 
D'yakonov-Perel' spin relaxation mechanism, which rep- 
resents a Markovian process, is the most important one 
in the wide range of t emp eratures and carrier concentra- 
tions. 01 00, Bill El EI E| At given sample parameters 
there are two possibilities to control the spin relaxation 
rate with this mechanism. First, one may apply an ex- 
ternal electric field in order to tune the Rashba constant 
0E1- Second, an external magnetic field can suppress 
the spin relaxation due to the cyclotron rotation, Larmor 
spin precession or their combination. [Til Hfij Although 
the main mechanisms of the magnetic field effects on the 
Markovian spin relaxation are understood, [lij j the pecu- 
liarities in the spin dynamics caused by the fascinating 
non Markovian transport phenomena in magnetic field, 
lavished attention only very recently [hR fisl flflj . and are 
scarcely known. Non-Markovian spin dynamics can arise 



for electrons interacting with nuclear spins [2(1 ] . hovewer, 
with a qualitatively different mechanism of the memory. 

In this Letter we point out and study a new non- 
Markovian effect in electron spin relaxation in the mag- 
netic field, namely, the appearance of the long-living spin 
polarization tail. We consider a degenerated 2DEG with 
SO coupling scattered by the short-range impurity cen- 
ters (or antidots) with the radius a and the concentration 
N in the presence of the magnetic field B directed per- 
pendicular to the structure plane. Assuming specular 
reflection by the scatterers, the kinetics of electrons is 
characterized by the total scattering length 1 = 1/ (2Na) 
and the transport mean free path l tr = 31 1 4. As shown 
in the pioneering work j2jj (see also Refs. I22I Eii . Eil and 
references therein) a fraction of electrons P C (B) moving 
with the velocity v does not experience scattering trav- 
elling along the closed circular orbits. In the Na 2 <C 1 
limit it is given by 



P c (B) = exp (-2ttR c /1), 



(1) 



where R c = v/lu c is the cyclotron radius, lo c = eB/mc 
is the cyclotron frequency, e is the elementary charge 
and m is the electron effective mass. Physically, Eq. |QJ 
gives the probability to find no scatterers in the area 
between circles of the radii R c ± a, see Fig. ^ As the 
electron finished the cyclotron revolution without a colli- 
sion it will stay on this orbit infinitely long, thus demon- 
strating a non Markovian behavior with an infinitely long 
memory. The fraction of these electrons increases with 
the magnetic field resulting in the classical magnetoresis- 
tance |23j. Here we discuss the role of the closed orbits 
in the spin relaxation. We show that the spin dynamics 
of the electrons on such trajectories is strictly periodical 
and the total spin falls exponentially down to the value 
determined by the fraction of circular orbits. Therefore, 
a tail in the spin polarization can be observed. Further, 
we develop a theory describing a decay of this tail caused 
by the scattering processes which transfer electrons be- 
tween closed and open "wandering" orbits. Finally, 
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wc discuss the role of the quantum effects. 
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FIG. 1: Electron's motion along the circling trajectory. Black 
circles - scattering centers. Solid line - initial "circling" tra- 
jectory. Dashed line - wandering trajectory after the inelastic 
scattering event with the scattering angle •&. Here a is the 
scatterer radius and b ~ N~ 1,/2 is the mean distance between 
scatterers. Figure shown not to scale. 

In order to obtain quantitative results wc consider a z || 
[001]-grownzince-blende quantum well (QW) and assume 
that the SO interaction is dominated by the Dresselhaus 
term Q 

^so(fe) = -^{c ■ J~Jfc), fife = ^lo[coatp k ,~sinif k ,0}. 

(2) 

where <r = (a Xl <j yi <j z ) is the Pauli matrix vector, Slo is 
the spin precession frequency at the Fermi level and (fk 
is the angle between electron wavevector k and x || [100] 
axis. 

We assume classical dynamics of the electrons and their 
spins and subdivide the electrons into two groups. First 
one corresponds to wandering electrons while the second 
one corresponds to circling electrons, their relative frac- 
tions arc 1 — P C (B) and P C (B) 1 respectively. The accurate 
description of the wandering electrons in the moderate 
magnetic fields can be done by using the kinetic equation 
and their spin decays exponentially with the rate |l4l | 
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t s {B) 1 + (w c t) 5 



(3) 



where r = Itr/vF- In the derivation of Eq. we as- 
sumed Hqt <C 1 ( 1 ) and neglected the Larmor effect 
on spin relaxation as it is usually small in GaAs-based 
QWs. The Larmor effect can be taken into account in a 
standard way. [Ir3 | 



1 This condition can be relaxed provided Qq <C u> c , M.M. Glazov, 
to be published 



The spin of the single electron with fc(i)-dcpcndent 
wavevector is described by equation 



— + s x n fc(t) 



o. 



(4) 



Eq. (0J can be readily solved assuming that k rotates 
in the QW plane with the frequency lo c . For the initial 
condition 2s(0) = (0, 0, 1), we obtain for z-component 



2s, (t) 
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V^o+^)- (5) 



The solutions for s XiV are rather cumbersome and not 
presented here. It is seen in Eq. JSJ that on the cir- 
cling trajectory 2s z oscillates around the mean value 
w c/( w c + fip ) ~± 1 ( m u ig u fields, u c » fio) with the 
frequency ^/fig + w 2 , and the amplitude of the oscilla- 
tions is fig/(w 2 + fig) — > at uj c ^> fi - Wc note that 
these oscillations are insensitive to the initial phase of the 
k rotation and thus all the spins of the circling electrons 
oscillate synchronously. 

Therefore, if in the 2DEG the transitions between cir- 
cling and wandering trajectories are strictly forbidden, 
at t 3> t s (B) (defined by Eq. ®) the total spin does not 
decay to zero as it was assumed previously, but oscillates 
around the non-zero value 



al = a P c (B)- 



+ fig' 



(6) 



where sq is the initially excited spin. Further we consider 
the case of the strong enough magnetic fields, uj c ^> fio 
where the oscillations of the circling electron spin are 
small and the spin relaxation is followed by the non- 
vanishing tail with the magnitude s* = P c (B)so. In order 
to have an insight into the tail formation we performed 
a Monte-Carlo simulation of the electron transport and 
spin dynamics with the results shown in Fig. [2J It is 
seen that the exponential decay of the velocity autocor- 
relation function K vv (t) = (v(i)v(0))/w 2 (shown in the 
inset) is followed by the oscillations due to the circling 
electrons. The magnitude of the oscillations corresponds 
to the value of P C (B). The polarization exhibits oscilla- 
tory tail in agreement with Eqs. (JSJ, ©. We mention 
that in nonquantizing field B = 0.5 T the tail has a 
considerable magnitude P C (B) > 0.1 already at low mo- 
bilities j( = 5x 10 4 cm 2 /Vs. 

Now we turn to the decay of the spin polarization tail. 
Any additional scattering process such as inelastic scat- 
tering or scattering by a random potential destroys the 
tail. With the relaxation time approximation the z spin 
component is described by the coupled kinetic equations: 



dsl s z [l-P c (B)]-s™P c (B) 

at ' ' Tj (b) 

ds? B ' z [l-P c (B)]-s?P c (B) 

at Tj(B) 



(7) 
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Here superscript c (w) corresponds to circling (wander- 
ing) electrons, respectively, and r s {B) and tj(B) is the 
spin relaxation time for wandering electrons and the in- 
termixing time between the closed and wandering tra- 
jectories. The first term in the numerator of the right 
hand side of the Eqs. 10 denotes the transfer of the elec- 
tron from the circling to the wandering trajectory and is 
proportional to the scattering intermixing rate and the 
fraction of the wandering trajectories 1 — P C (B) and the 
second term refers to the inverse process. In deriving 
Eqs. Q we have neglected small oscillations of s z and 
took into account D'yakonov-Perel' spin relaxation of the 
wandering electrons (second term of the second line in 
JjJ). The eigenrates of the system r-f^-B) are 

where T js = ^(B) + t- 1 (B). Thus, the total spin of 
electrons s z (t) relaxes according to the two exponential 
law: s z (t) — aexp (—t/T\[B)) + /3cxp {—t/rz^B)) where 
coefficients a and (3 depend on the initial condition and 
B. Since t 2 (B) > Ti(B) long-time s z (t) asymptotics is 
determined by T2(B) which can be identified with the tail 
relaxation time. 

For example, in the limit of weak magnetic fields 
P C (B) <C 1 and if t s (B) <C Tj(B) the relaxation times 
are T\{B) = t s (B) and t 2 (B) = Tj(B) and the total 
electron spin rapidly decays to the tail value P c (B)sq 
and then slowly relaxes during the time Tj{B) (as almost 
every " additional" scattering transfers electron from the 
circling to wandering trajectory where the spin is quickly 
lost). 

When the magnetic field becomes stronger with ui c t 3> 
1, t s (B) = r s (0) (u! c t) 2 increases quadratically with the 
field and the fraction P C (B) = 1 — 3tt/2lu c t tends to unity. 
Two limiting regimes as a function of B are possible. In 
weak fields, where the intermixing Tj is still much longer 



than t s (B), the spin tail relaxation time is field inde- 
pendent. With the increase of B, the intermixing time 
becomes smaller than t s (B), and, therefore, the tail dis- 
appears. In such a case only a small fraction (~ lo c t) of 
electrons is loosing at a given moment their spin, how- 
ever circling and wandering electrons are intermixed on 
a short timescale. Therefore the tail is not observed but 
the spin relaxation time is: 

being proportional to B 3 (contrary to the standard re- 
sult B 2 |l4|P- Here t 2 (B) does not depend on the inelastic 
scattering time despite the presence of the strong inelas- 
tic scattering which is necessary here only to intermix the 
trajectories sufficiently quickly. These results are shown 
in Fig. [3]a, illustrating the evolution of the spin relaxation 
process s z (t) and the rate with the change in the number 
of wandering trajectories governed by the independent 
ratio I / Re- 
in sufficiently strong fields, the number of trajectories 
colliding several times with the same impurity, increases, 
and at a critical field such as 4irN(a + R c ) 2 m 4.4 and 
(R c /l) 2 = 4ANa 2 /n < 1 HI] all the paths become con- 
fined. For the parameters considered in Fig. (2), this cor- 
responds to B > 0.8 T. In this case the spin relaxation 
is strongly suppressed j^. The detailed analysis of this 
regime is beyond the scope of this Letter. 

Below, we consider two kinds of the microscopic pro- 
cesses which contribute to the intermixing time Tj(B) 
between closed and wandering trajectories. Namely, we 
discuss the destruction of the tail due to: (i) presence of 
the smooth potential, (ii) inelastic processes. 

(i) In the real QWs the typical disorder can be con- 
sidered as the superposition of the short-range scattering 
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FIG. 3: (a) Electron spin s z time dependence calculated from Eqs. J7J for the different values of the magnetic field: l/R c = 1, 2, 5 
from the bottom to the top, Tj = 10r s (5 = 0). (b) The scattering probability found by the Monte-Carlo simulation for (solid 
lines) and values of P C (B) (dashed). The parameters used in the calculations are shown above each curve. Other parameters 
are same as in Fig. (2). The accuracy of the Monte-Carlo computation is shown by the error bars. 



centers ( 2 ) and the long-range smooth potential of the 
remote dopants and interface roughness ( 3 ). The correla- 
tion radius of this potential d is the distance between the 
dopant layer and the QW. [2(| ^ We will be interested 
in the weak smooth disorder where the transport scatter- 
ing time by long-range potential 77 satisfies the condition 
Titu c S> 1, and, therefore, the circling trajectories remain 
at least approximately intact. The displacement of the 
orbit center in the smooth potential provided R c d 
is Sr ~ k/iuj^f/ 2 ~ ^/Rj/d(W/E F ), where h = v F 
r\ is the transport mean free path there, W is the po- 
tential amplitude. [2^] The transfer between trajectories 
is governed by the ratio Sr/a. The rich transport [2(| 
and spin dynamics of the electrons in the system with 
two types of disorder are beyond the scope of this Letter. 
Here we consider the simplest situation where Sr <^ a 
and thus Sr <C d. Therefore, electron drifts in the poten- 
tial gradient averaged over its path and sweeps during n 
revolutions an area S ~ 27r R c n5r. It will be transferred 
from the circling trajectory to the wandering as soon as 
it meets a short-range center in this area, i.e. SN ~ 1. 
So, the intermixing time can be estimated as 

T ? {B) ~ m& « s3/2 ' 6r/a<<1 > (10) 

where T c = 2tt/lu c is the cyclotron period and it is as- 
sumed that alt r /(dR c ) <C 1. For the purpose of the 
present paper we note that the electron drift in the 



2 These cent ers can be introduced artificially as antidotes, see e.g. 

03 nam 

3 The monolayer fluctuations of the QW interfaces induce the po- 
tential 5V ~ 1 . . . 10 mcV with the lateral size 50 . . . 100 A. Thus, 
the scattering length on the interface roughnesses a is of the order 
of 10 4 A 



smooth potential leads to the spin relaxation even for the 
circular trajectories, therefore the applicability of the Eq. 
(tTU)l together with the system (jTJ) is guaranteed provided 
if the intermixing time is much longer than the spin re- 
laxation time of the circling electrons. 

(ii) The inelastic processes (electron-electron or 
electron-phonon scattering) change randomly electron 
wavevector and energy and thus transfer electrons be- 
tween circling and wandering trajectories. In the low 
temperature regime where the Fermi surface is well de- 
fined, Ep 3> fcsT, where Ep is the Fermi energy and 
ksT is the temperature measured in the energy units all 
the scattering processes are statistically quasi-elastic , 
i.e. the transferred energy ~ ksT is much smaller than 
electron kinetic energy Ep. The crucial question is to 
determine the probability that the scattering at the an- 
gle $ (Fig. HJ is not accompanied by the electron tran- 
sition from the circling to wandering orbit. From Fig.^ 
it is clear that the displacement of the circular orbit by 
the distance of the order of the impurity separation b 
will lead to the electron collision with some scattering 
center thus the probability should have a peak in the 
vicinity of $ = with the angular width of the order 
of b/R c - \/{N 1 / 2 R c ). At lar ger angles the probabil- 
ity saturates at the value close to P C {B) representing the 
fraction of the circling electrons. This prediction is cor- 
roborated by the Monte-Carlo calculation, see Fig. 03. 

In a degenerate 2D electron gas with carrier concen- 
tration ranging in 10 11 . . . 10 12 cm~ 2 the Coulomb inter- 
action is screened and the electron-electron scattering is 
almost isotropic and its temperature dependency is de- 
termined by the Pauli-exclusion principle only. At not 
too low temperatures fcgT > hc s \J (fc 2 ) where c s is the 
speed of sound, and y (^ 2 ) is the root mean square of the 
electron wavevector in the growth direction, the electron- 
acoustic phonon scattering is dominated by deformation- 
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potential coupling and is isotropic as well. Thus we fur- 
ther assume that the inelastic scattering is isotropic and 
can be described in the framework of the relaxation time 
Tj which can be estimated as for electron-electron and 
electron-phonon scattering, respectively [3(il.l3lj |: 

ee_ r J_f^F_\ 2 ph_, PC 2 s h 3 

~^ e E F \k B T) ' ~~ D 2 mk B T«/(k%}' 

(11) 

Here ( e — A/[tt\ii(Ef /ksT)] for Fermi level electrons, 
D is the deformation potential constant, p is the ma- 
terial density, and ( p h ~ 1 is the numerical coefficient. 
At lower temperatures the electron-phonon scattering is 
partially suppressed due to Pauli-cxclusion principle and 
by the phonon-freeze-out and the coefficient Cph becomes 
strongly temperature dependent. Its exact calculation 
represents an extremely difficult numerical tasks but the 
simple estimates show that in the temperature range 
T < 20 K and for electron concentration N e ~ 10 11 
cm" 2 the electron-electron scattering is dominant. We 
note that the fraction of the circling electrons P C (B) Eq. 
(|T|) has a purely geometric character, thus the tail is in- 
sensitive to the thermal broadening of the electron dis- 
tribution provided inelastic processes are weak enough. 

Finally, we would like to comment on the quantum ef- 
fects neglected above. First, the quantization of electron 
orbits leads to the oscillatory independence of the spin 
relaxation rate (0) [13 ■ This effect is small for high 
Landau levels (Ep/huj c 3> 1) and vanishes as tempera- 
ture increases. Second, the quantization of the electron 
orbits modifies the relaxation rates of inelastic scattering 
in Eq. i|ll|) . However, for semiclassical Landau levels 
this modification is negligible [i^. Third, the magnetic 
field leads to the spin polarization cither due to Zeeman 
effect or as a result of the SO coupling and Landau quan- 
tization 0. The former contribution is proportional to 
gfjLsB I Ep where g is the electron g- factor and [is is Bohr 
magnetron. For the typical GaAs structures g is small, 
besides it can be adjusted to zero with the structure pa- 
rameters @- The contribution \Vl\ = \{4 , l,±Wz\4>l,±)\ 
of the latter effect can be evaluated for Lth Landau level 
E L> ± = hio c (L± ^7 2 L+ 1/4) with 7 = hn /(2y/E F huj c ) 

as = (2-\/7 2 L + 1/4) 1 and is negligible for high 
Landau leves with 7 2 L> 1. 

In conclusion, we have theoretically analyzed the spin 
dynamics of the 2D electron gas scattered by short-range 
defects in the classically strong magnetic fields. We have 
shown that the tail in the spin polarization appears as a 
result of the collisionlcss circular motion of the fraction of 
electrons. It was demonstrated that the tail decays due 
to the presence of the additional weak long-range poten- 
tial or due to inelastic scattering. In strong magnetic 
fields the long-time spin relaxation is slowed down as a 
cube of the field. The predicted phenomena open new 
perspectives in the spin dynamics control in the 2DEG. 
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